Sequential Tunneling through Molecular Spin Rings 
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We consider electrical transport through molecules with Heisenberg-coupled spins arranged in 
a ring structure in the presence of an easy-axis anisotropy. The molecules are coupled to two 
metallic leads and a gate. In the charged state of the ring, a Zener double-exchange mechanism 
links transport properties to the underlying spin structure. This leads to a remarkable contact-site 
dependence of the current, which for an antiferromagnetic coupling of the spins can lead to a total 
suppression of the zero-bias conductance when the molecule is contacted at adjacent sites. 

PACS numbers: 85.65.+h, 85.75.-d, 75.10.Jm, 73.63.-b 



During the last decade, we have witnessed a tremen- 
dous progress in the experimental methods for contact- 
ing single molecules and measuring the electrical current 
through them [1-5]. A break-through has been achieved 
with the electromigration-j unction technique [2], which 
allowed one to add a back-gate to the molecule, thereby 
bringing the whole field of Coulomb-blockade physics to 
molecular systems. Recently, this technique has led to 
the first transport measurements through single magnetic 
molecules: The current through the different Mni2 clus- 
ters [6] has been measured and magnetic excited states 
have been identified [3, 4]. At the same time, various 
theoretical models for the description of these measure- 
ments have been put forward, all based on the effective 
Hamiltonians routinely used for the description of ther- 
modynamic properties of the molecular magnet [3, 4, 7- 
9]. Doing so, an immediate complication encountered is 
that for the description of the sequential tunneling pro- 
cesses observed in the experiment, at least one charged 
state of the molecule has to be included in the model, 
as well. As plausible models for this state, the effective 
Hamiltonian with renormalized parameters [3, 7] or mod- 
ified by a single additional orbital [8] have been employed. 

Here, we take a microscopic approach in terms of the 
constituent spins of the molecular magnet. Whereas this 
limits us to smaller systems than the ones considered ex- 
perimentally, it enables us to include an excess charge. 
This will shed some light on the role of the orbital de- 
gree of freedom of this excess charge (say an electron) 
not contained in an effective-Hamiltonian model. Fur- 
ther, we are able to consider on the same footing trans- 
port through antiferromagnetically (AF) coupled rings 
like molecular ferric wheels [10], which has not been stud- 
ied so far. 

Evidently, an excess electron will reduce one of the 
positively charged ions forming the localized spins. The 
arbitrary choice of this ion leads to an orbital degen- 
eracy which will, however, be lifted by a hopping pro- 
cess between the different sites. At the same time, this 
hopping allows the electron to travel from one side of 
the molecule to the other and a current to flow across 
the molecule. Since according to Hund's rules, the cou- 
pling of the excess electron to the localized spin will favor 



one direction of the excess electron's spin, and assuming 
that hopping leaves the spin unchanged [17], one obtains 
an effective ferromagnetic (FM) type of coupling known 
as Zener double exchange [11]. Correspondingly, the or- 
bital degree of freedom of the electron is linked with the 
spin structure of the ring. In particular for AF rings 
this should lead to interesting transport effects: Due to 
the Neel-ordered structure of the AF ground state, we 
even expect — and indeed confirm — a strong even-odd- 
site-dependence of the current mediated by ground-state 
transitions not present in the FM case. Another in- 
triguing property of AF rings is that quantum tunneling 
of magnetization (QTM) leads to a tunnel-split ground- 
state doublet with a characteristic oscillating dependence 
on the magnetic field [12]. We will demonstrate how this 
behavior can be extracted in a transport measurement. 

Model — The physical picture described in the previous 
paragraph is captured by the Hamiltonian 
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Here, the first term describes the Heisenberg coupling 
between the localized spins Si with spin quantum num- 
ber s, where we require periodic boundary conditions 
Sat+i = Si to describe a ring with total spin S = ^ 
The second term represents the Hund's rule coupling of 
the excess spin when it is localized at the site i to the 
corresponding spin S^. The second-quantized spin oper- 
ators Si are defined in terms of the creation (annihilation) 
operators d\ a for an electron of spin a at site i by 

Si = (1/2) Y^aa' T ™' d l,a d hv' witn T being the vector of 
the three Pauli matrices. The spin of the excess electron 
is then given by s = X^ s *- The easy- axis anisotropy 
of strength k z and the Zeeman energy are included by 
the third and fourth term, respectively. The fifth term 
describes the hopping of the excess electron around the 
ring, where again periodic boundary conditions are as- 
sumed. The orbital energy eo of the excess electron can 
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be tuned by applying an external gate voltage V g and is 
proportional to the occupation number n = ^ d\ G di^ . 
Finally, we restrict ourselves to a single excess electron 
by letting the charging energy U — > oo. 

The coupling of the molecule to the metallic 
leads is described by the tunnel Hamiltonian Ht = 
X^kcr cJ ko .d^a + h.c. Here, the sum runs over the two 
leads i — L, R, which we assumed to be each coupled via 
a spin-independent coupling matrix elements to a sin- 
gle ring site ig. The leads at an electro-chemical potential 
fii are modeled as a set of independent quasiparticles with 
wave-vector k and spin a, H L = ^ ih(J Q k cJ ka Q ka , dis- 
tributed according to the Fermi distribution f(e — = 
{1 + exp[(e — fii) //cT]} _1 . Assuming that the exter- 
nally applied bias voltage Vb drops completely and sym- 
metrically along both molecule- lead contacts, we have 
Ml,r = /i ± eVb/2, where \i is the Fermi energy. 

Magnetic order — The eigenstates \a) of the unper- 
turbed ring Hamiltonian (1) can be separately considered 
for the cases of an uncharged (n = 0) and charged (n = 1) 
ring. Moreover, the discrete translational symmetry and 
the symmetry of reflection at an arbitrary site, say i = 1, 
yield corresponding good quantum numbers cos(q) with 
— 7r < q < 7r and Nq = 0mod27r, and parity P = ±1, 
respectively. In the absence of a magnetic field with x- 
and ^/-components, another good quantum number M is 
given by the z-component of the total angular momen- 
tum S + s. The transition matrix elements (a f \d\ a \a) of 
the creation operator fulfil the selection rules n' — n = 1 
and M r — M = a. In the special case N = 4 considered 
below, we find P' = P for odd site numbers i = 1, 3 and 
exp(2V)P / = exp(2ig)P for even site numbers i = 2,4. 

We now discuss the magnetic order of the lowest-lying 
states for the two opposite cases of a FM and AF intra- 
ring coupling J. In both cases, we assume that the ions 
forming the localized spins have at least half-filled shells 
and thus consider according to Hund's rules a strong AF 
coupling Jh <C — \J\. This describes, for instance, ferric 
wheels containing Fe 3+ ions. We will exemplify the gen- 
eral discussion by the special case of a ring consisting of 
N = 4 spins s = 1 (cf. Fig. 1). The exchange coupling 
strength \J\ serves as energy unit (typical values for ferric 
wheels range from 20 to 30/cbK [10]). The anisotropy is 
set to k z = 0.3 1 J\ and the Hund's rule coupling strength 
is Jh = — 100 1 J|. An estimate for the hopping matrix 
element t is difficult to make. Assuming that it is of the 
same order of magnitude as the one leading to exchange 
in the uncharged ring, a simple argument based on a 
Hubbard model yields the lower bound J /At « t/U' <C 1, 
where U 1 is the on-site charging energy. Expecting the 
excess electron to be rather localized, we assume t = 5|J|. 

For a FM intra-ring coupling, the lowest-lying states of 
both the neutral and charged ring consist of multiplets 
of total spin S — Ns and S — Ns — 1/2, respectively 
[cf. Fig. 1(a)]. These multiplets are split by the easy-axis 
anisotropy into doublets with same absolute value of the 
magnetic quantum number \M\. The additional FM-type 
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FIG. 1: (Color online) Low-energy states of the neutral (left, 
n = 0) and charged (right, n = 1) ring both for (a) FM 
and (b) AF coupling. The parameters are given in the main 
text. The levels are labeled by their magnetic quantum num- 
ber M and their parity P. Dashed (dotted) lines indicate 
tunnel transitions allowed when the ring is contacted at the 
odd (even) numbered sites. Transitions allowed for all contact 
sites are drawn as dash-dotted lines. 



interaction due to the Zener double exchange influences 
the magnetic order only weakly. 

This is in stark contrast to the case of an AF intra- 
ring coupling J < [cf. Fig. 1(b)], which, for the charged 
ring, competes with the double-exchange contribution. 
In order to understand the resulting magnetic ordering 
in this case, we first consider the uncharged ring. The 
classical ordering is described by the two Neel states 
\t) = |s, -5,5, -s, ...) and = 1-5,5,-5,8,...), 
where the vectors have been written in the {S z ^} eigen- 
basis. In the QTM reg ime de fined by the scaled tun- 
nel action So/h = Ns^/2k z / J being much larger than 
unity [12], the leading contribution to the two lowest- 
lying quantum-mechanical eigenstates of the uncharged 
ring is then given by the symmetric and antisymmetric 
superposition of the two Neel states, i.e., (|ff) ± \§))/y/2. 

Similarly, in a purely classical picture and without hop- 
ping contribution, the charged ring has a 2N-io\d de- 
generate ground state which for an AF on-site coupling 
Jh < consists of the states |ff) \i, j) and |J|) \i, j) for odd i 
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as well as |) and |J|)|i,t) f° r even Here, |i,<r) de- 
scribes the state of the excess electron being localized at 
site i with spin a. The hopping leads to an hybridization 
of the states of same Neel ordering and excess electron 
spin. Furthermore, QTM couples states with opposite 
Neel ordering of the localized spin lattice. 

We find for the special case with N = 4 and 5 = 1 
[cf. Fig. 1(b)] that the ground and first excited state of 
the uncharged ring are indeed dominated by the super- 
position of the Neel-ordered states described above. The 
ground state of the charged ring turns out to be four- 
fold degenerate, discerning into two degenerate Kramers 
doublets with M = ±1/2. Focusing on the states with 
M = 1/2, we find that in the parity (P) eigenbasis, 
the main contribution to the eigenstates stems from the 
states |fr)(|2,|) - |4, |)) with positive parity P = 1 
and |J|)(|1,|) — |3,|)) with negative parity P = —1. 
These states do not exhibit a QTM-induced tunnel split- 
ting [18]. The situation is different for the two next 
excited Kramers doublets. The leading contribution to 
these states is given by the symmetric and antisymmet- 
ric combination of the two states j) + |3,|)) and 
|J|)(|2, t) + |4,|)) and we indeed observe a finite tunnel- 
splitting due to QTM. 

Block- Redfield approach — We calculate the sequential- 
tunneling current through the ring using a Bloch- Redfield 
approach, extending previous results [13]. In contrast 
to a Pauli master-equation description, this includes off- 
diagonal elements (coherences) of the reduced density- 
matrix g a (3 = Tr^, (a\g\f3) after tracing out the leads and 
is also valid for level spacing less than the level broaden- 
ing induced by the leads. We find 
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Here, we have introduced the frequencies u a p = 
(E a — Ep)/h and the transition rates Wp,p aa , = 

W^ aa , + Wp7p aa , due to tunneling across a molecule- 
lead contact £ in terms of the rates Wptp aa t = 
T £ (E a -E af ) f(E a -E a ,-»i) E CT 

and W P 7 P**> = Wa> ~ E a ) [1 - f(E a , 

Vi)] ^2a(^ f \ d lia\^)( a \ d U^W) for tunneling of an elec- 
tron on and off the molecule, where T^(e) = 
(2ty/H) ^ k \T^\ 2 5(e — e^k) is the spectral density of the 
coupling. Finally, the current across contact £ is given 
by/, = e Re£ a ^«^a-< + 



The number of rates Wp f/3aa , 
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scales with the fourth 
power of the dimension of the molecular Hilbert space. 
For the numerical solution of Eq. (2), we thus cannot take 
into account all coherences. However, for the station- 
ary solution, coherences between states with an energy 
difference much larger than the molecule-lead coupling 
strength are negligible. 

Differential conductance — For our numerical calcula- 
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FIG. 2: (Color online) Differential conductance (color-coded) 
as a function of gate and bias voltage at zero magnetic field. 
The temperature is k^T — 0.05 \J\ and the other parameters 
are as given in the main text. The upper two panels (a) and 
(b) show the conductance for an FM coupling J > 0, the 
lower panels (c) and (d) depict the AF situation J < 0. The 
left (right) panels show the conductance when the molecule is 
contacted at opposite (adjacent) sites = 1 and ir = 2 (zl = 
1 and in = 3) (cf. inset schematic drawings). The dashed 
lines indicate the position of the ground-state transitions. 



tions, we use the model parameters described above and 
assume a contact- and energy-independent molecule-lead 
coupling ftr = hT£ = 0.01 1 J|. Figure 2 shows the differ- 
ential conductance G = dl/dVb as a function of bias and 
gate voltage at zero magnetic field for both the FM and 
the AF case and for two molecule- lead coupling variants: 
contacts at opposite and adjacent, respectively, sides of 
the ring. For a FM coupling J > 0, the behavior in 
the low-bias- voltage regime is identical in both cases [see 
Figs. 2(a) and (b)]: the ground-state transitions domi- 
nate the current and one excited state can be identified. 
Away from the charge-degeneracy point, one ground- 
state transition is suppressed for adjacent contacts, for 
reasons similar to the AF case discussed next [14]. 

In the case of an AF intra-ring coupling J < 0, we find 
a more striking difference between the two contacting 
situations. For adjacent contacts [see Fig. 2(d)], the low- 
bias conductance is completely suppressed, which does not 
occur for contacts at opposite sides [see Fig. 2(c)]. In or- 
der to understand this remarkable difference, we consider 
the magnetic order discussed above. There, we identified 
two Kramers doublets of different parity in the n = 1 
ground state. Furthermore we found the contact-site de- 
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FIG. 3: (Color online) (a) Differential conductance through 
an AF coupled ring contacted at opposite sides as a function of 
magnetic field in x-direction and bias voltage. The gate volt- 
age has been chosen as indicated by the arrow in Fig. 2(d). 
The temperature is k&T = 0.01 \J\ and the other parameters 
are given in the main text. The solid (dashed) lines indicate 
the electron (hole) transition labeled in panel (b). The tun- 
nel splitting A extracted from these transitions is shown as 
white dotted line, (b) Level structure of the ring and tunnel 
transitions for finite B x as indicated by the arrow in panel (a). 

pendent selection rules depicted in Fig. 1. In particular, 
when contacting at adjacent sites, one contact only al- 
lows transitions to the even parity states and the other 
one only to the ones with odd parity. Thus, there is 
no transport path connecting both contacts and the cur- 
rent is blocked completely until transitions via excited 
states become possible at higher bias- volt ages. On the 
other hand, for contacts at opposite sides, only one par- 
ity state becomes populated and a current can flow. It 
displays a rather complex excited-state spectrum com- 
pared to the FM case. For a better understanding of this 
regime, we focus on a fixed gate voltage as indicated by 
the arrow in Fig. 2(c) and consider in Fig. 3(a) the con- 
ductance as a function of an in-plane magnetic field. This 



dependence is particularly interesting for AF rings where 
the QTM-induced tunnel splitting A between ground and 
first excited state of the neutral ring [see Fig. 3(b)] os- 
cillates as a function of the magnetic field [12]. In or- 
der to extract this splitting from the conductance behav- 
ior, we identify all transitions between the lowest-lying 
states of the neutral and charged molecule (see labeled 
lines in Fig. 3). Denoting the magnetic-field dependent 
bias voltages at the corresponding resonances by Vh,k 
(k = 1,...,4), the tunnel splitting can be obtained as 
A = e(Vb,i + Vb, 3 )/2 = e(V ha + V hA )/2. 

We finally remark that a similar contact dependence 
of the current, albeit not a complete suppression of the 
zero-bias conductance, has been predicted for benzene 
rings [15]. There, it was also noted that when excited 
states are populated, a non-tunneling decay may lead to a 
population of the odd-parity state such that for contacts 
at opposite sides one may find a current collapse at higher 
bias voltages. A similar scenario can be envisaged here 
with the additional possibility to lift the blockade again 
by applying a large-enough magnetic field B x [14]. Then, 
the current should display an interesting hysteric behav- 
ior possibly useful for information storage purposes. 

To conclude, we have presented a model for a molecu- 
lar spin ring in the presence of a single excess electron. It 
contains a hopping term which both provides a transport 
path and results in a Zener double-exchange coupling be- 
tween the localized spins. For AF coupled rings of ions 
with least half- filled shells, the interplay of these two ef- 
fects leads to a strong contact-dependence of the sequen- 
tial tunneling current through the molecule reflecting the 
Neel-ordered ground states. 
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